Atomic quantum simulator for lattice gauge theories and ring exchange models 
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We present the design of a ring exchange interaction in cold atomic gases subjected to an optical 
lattice using well understood tools for manipulating and controlling such gases. The strength of 
this interaction can be tuned independently and describes the correlated hopping of two bosons. 
We discuss a setup where this coupling term may allows for the realization and observation of 
exotic quantum phases, including a deconfined insulator described by the Coulomb phase of a three- 
dimensional U(l) lattice gauge theory. 



Loading cold atomic gases into optical lattices allows 
for the realization of bosonic and fcrmionic Hubbard 
models, and offers the possibility for the experimental 
study of strongly correlated systems within a highly tun- 
able environment. Starting from the prediction of a su- 
perfluid to Mott-insulator phase transition in bosonic 
atomic gases and the subsequent observation of the 
Mott insulating phase [3- 0- [~3] - many new tools for ma- 
nipulating and controlling quantum gases have been de- 
veloped |5j]. In this Letter, we combine these tools in 
order to drive the atomic gas with an additional ring 
exchange interaction. We identify a promising system 
where this coupling may allow for the realization and ob- 
servation of an exotic quantum insulator 0, described 
by the Coulomb phase of a three-dimensional U(l) lat- 
tice gauge theory ||; in quantum magnetism this phase 
is known as a U(l) spin liquid. 

Recently, studies of boson models with large ring ex- 
change have yielded significant progress in the search 
for microscopic Hamiltonians exhibiting exotic phases 
0- fiol HH - This search has been the focus of much 
effort in two-dimensional systems, due to potential rel- 
evance for high-T c cuprates. Some ring exchange mod- 
els exhibit a local conservation law, and can be mapped 
onto lattice gauge theories and often also quantum dimcr 
models (QDM) [l^. A number of such models in two 
and three dimensions have been shown to possess de- 
confined insulating ground states Many of the 
three-dimensional models, including those of Refs. 0, 0, 
were shown to possess a U(l) deconfined phase, which 
supports gapped half-boson excitations, gapped "mag- 
netic monopole" topological defects, and a linearly dis- 
persing photon mode with two polarizations. The low- 
energy theory is standard quantum electrodynamics with 
massive electrically and magnetically charged scalar par- 
ticles. Models of bosons on the square lattice with 
large ring exchange are in a different class from those 
above, and also exhibit interesting physics. Such mod- 
els can support an "exciton Bose liquid" phase, a two- 
dimensional analog of a Luttinger liquid as well 



as nontrivial valence-bond solid (VBS) insulating states 
[lfij . Such states can undergo a direct quantum phase 
transition to the superfluid 16]. 

Despite much recent theoretical progress, clear experi- 
mental evidence for the existence of exotic phases is still 
missing. Furthermore, relatively few theoretical tech- 
niques exist to study such strongly correlated systems; 
perhaps the most powerful to date is quantum Monte 
Carlo simulation, but the class of models that can be 
productively studied is severely limited by the notorious 
sign problem. Atomic gases offer an alternative approach 
through the design of quantum simulators, where a mi- 
croscopic Hamiltonian is implemented in a quantum gas 
and its phase diagram is studied experimentally via con- 
trolling the strength of the interaction terms. 

In this Letter, we present the design of a ring exchange 
interaction for bosonic gases subjected to an optical lat- 
tice. Such an atomic lattice gas is well described by the 
Bose-Hubbard model Q 



H BH = -jY,blb J + ^J2 b ^ b 



(1) 



(id) 



where U denotes the on-site repulsion and J is the hop- 
ping energy with denoting summation over nearest- 
neighbor sites. The additional ring exchange interaction 
involves four lattice sites forming a square plaquette, and 
is driven by a resonant coupling of the bosons via a Ra- 
man transition to a "molecular" two-particle state jjj, 
see Fig. This state is subjected to an independent op- 
tical lattice with its lattice sites at the plaquette centers. 
The symmetry of the molecular state strongly influences 
the coupling; we are interested in a d-wave symmetry of 
the molecule, which can be carried either by the relative 
coordinate or the center of mass motion. Then, the cou- 
pling to the molecular state (created by nv a ) takes the 
form 

H M = ^2vmlm n +g^2 [ m a i b ih - b 2 h) +R.c] . (2) 

□ □ 

The summation runs over all plaquettes □. The single- 
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particle states created by b\ are called bosonic or atomic 
states to distinguish them from the "molecular" two- 
particle states. Depending on the setup, the atomic 
states reside either in the corners or on the edges of each 
plaquette (see Fig. Q), and are numbered counterclock- 
wise. The energy v corresponds to the detuning from 
resonance, while g is the coupling strength determined 
by the Rabi frequency of the Raman transition. While 
the Hamiltonian (J2J) is interesting in its own right, the 
connection to ring exchange is apparent upon integrat- 
ing out the molecular field perturbatively in g/v, which 
leads to the effective Hamiltonian 



■nin 3 -n 2 n i 



(3) 



with K = g 2 jv. Note that the structure of the coupling 
in Eq. (0 also produces a next-nearest-neighbor inter- 
action. The bosonic system turns metastable for large 
negative detuning. However, the decay time easily ex- 
ceeds typical experimental time scales of atomic gases. 
Then, the perturbative expansion is again valid and al- 
lows for the realization of a system with negative ring 
exchange interaction. 




FIG. 1: (a) Two-dimensional setup: the bosons (black dots) 
are on the square lattice with the molecules (gray dots) in 
the center of each plaquette. (b) Three-dimensional setup: 
the bosons (black dots) are on the links of the cubic lattice. 
Within each face there are four bosonic sites, which establish 
a plaquette (dashed square). The molecules (gray dots) are 
in the center of each plaquette. 



The design of the ring exchange Hamiltonian H = 
H Bn + H M combines standard tools for manipulating and 
controlling atomic gases: First, we are interested in a 
cold atomic gas with two internal states, which are sub- 
jected to independent optical lattices; the internal states 
can be coupled via a Raman transition. Such a setup 
has been realized recently using spin-dependent optical 
lattices 0, 0] . An alternative approach is the trapping 
of alkaline earth metals, e.g., 88 Sr, where the first ex- 
cited state 3 Pi exhibits a long life-time with a different 
polarizability than the lowest energy state 1 5o 0]. Sec- 
ond, combining different laser configurations allows for 
the design of complex lattice structures beyond the stan- 
dard cubic lattice |20| , and third, the interaction between 



the atoms can be tuned by magnetic or optical Feshbach 
resonances |5|. 

We introduce the notation ^> a (x) and ip b (x) for the 
field operators describing the two internal states. The 
microscopic Hamiltonians then takes the form (a = a, b) 



dx 



with g a = Aith 2 a a /m the interaction strength for scat- 
tering lengths a a . The e a are the homogeneous en- 
ergy shifts between the internal states. The potential 
V a (x) accounts for an optical lattice driven by lasers with 
wavevector k = 2tt/X, with the strength v a in units of the 
recoil energy E r — h 2 k 2 /2m. The two internal states are 
coupled via a Raman transition. Transforming away the 
optical frequencies within a rotating frame, the coupling 
takes the form 



(4) 



with Q the Rabi frequency of the transition. 

We focus first on the two-dimensional setup shown 
in Fig. ^i. Confinement to two dimensions is achieved 
by a strong transverse optical lattice, which quenches 
hopping between different planes. The remaining op- 
tical lattice provides the square lattice structure for 
the atomic state -06; an d takes the form V h /E r = 
v b [cos 2 (kx/2) + cos 2 (fcy/2)] . For v b > 1, the mapping 
to the Bose Hubbard model is well justified. The optical 
lattice for the second internal state i/> a , which is localized 
at the plaquette centers, takes the form 



v a ([coskx— cos/cj/] 2 + sin 2 (fcx/2)+sin 2 (fcy/2) j 



The first term is obtained by interference between stand- 
ing laser waves along the x- and y-directions, while the 
other terms represent a standard square lattice. The dif- 
ferent lattice spacing of the two contributions is easily 
achieved by a finite angle 2-7r/3 between the interfering 
beams. We are interested in a strong optical lattice V a , 
where tunneling between different wells is quenched, and 
focus on the energy states within a single well. Then the 
structure of V a produces strong shifts of the energy states 
compared to those obtained within a harmonic approxi- 
mation. The states are characterized by the irreducible 
representations of the symmetry group Ci V {i.e. the 
point group of the square lattice); the low-energy states 
and corresponding representations derived within a band 
structure calculation are shown in Fig. [21 for v a = 30. 
The state with energy e; in each plaquette □ is created 
by the bosonic operator a\ a with I = 0, ±1,2. Of par- 
ticular interest is the state |a2,n) corresponding to the 
representation B 2 , which transforms under d v like the 
polynomial xy {i.e. d xy symmetry). In contrast to the 
harmonic approximation, this state is non-degenerate. 
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For weak interactions, the Hamiltonian for the bosonic 
field -0a reduces to 



H„ = 



l, a 



via\ a i 



■i,n 



E 



Ui,v 
2 



l La a i'a a l,a a i',o 



(5) 



with Uiji < ei the interaction energy, and vi = e/ — fko 
the energies of the excitations within the rotating frame 
(lu is the frequency of the Raman transition). The cou- 
pling driven by the Raman transition takes the form 



with £1 the Rabi frc- 



quency. Due to the square symmetry, each operator oj t n 
couples to a special structure of surrounding bosons. The 
operators for which the coupling becomes diagonal trans- 
form irreducibly under C^ v and are do,n ~ 61+62+63+64, 
d±\,u ~ 6i±i6 2 — 63=^64, andd 2 ,n ~ 61—62+63 — 64. The 
wave function overlaps wi derive from the shape of the 
localized wave functions. For typical parameters (v h ~ 6 
and i> a ~ 30) we obtain wi w 0.1. 
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FIG. 2: (a) Single-particle energies e; of the states a; n |0) 
for v & — 30. (b) Energy levels of the two-particle states with 
symmetry B2 and E. The frequency lo = (eo+£2+£/o2+^)/2ft 
of the Raman transition is chosen near resonance with the 
molecular state aJa^O), and far-detuned to the single-particle 
excitations. 



We are interested in a setup with "molecular" two- 
particle states \m a ) = Wp|0) = aJc^onlO) (with d xy 
symmetry B2) resonantly coupled to the bosonic states 
6t £>t j 0) with detuning v. The energy of this state is e M = 
£0 + £2 + Uq2 ■ As the formation of the molecule involves 
the virtual creation of a single particle excitation at n |0), 
the frequency of the Raman transition is determined by 
lu = (e M + u)/2h (Fig. EJ). We require that the single- 
particle excitations a^~ n |0) are far-detuned {i.e. wiMl < 
\vi\), and only the molecular states are resonantly coupled 
with \v\ -C As the separation between the energy 
levels e; is large (~ 2-Jv^E r ), this condition can easily 
be satisfied for suitable interaction energy Uo2- Further- 
more, it may be of interest to tune the interaction en- 
ergy via a Feshbach resonance to an optimal value. For 
such strong interactions the molecular operator becomes 



mt =cia6,p a 2,p+ c 2( a i,nai,n +a -i,n a -i,n)" 1 ' where 

the ellipsis denote admixture of higher energy states re- 
specting the d xy symmetry. The parameters c\p and the 
energy e M have to be determined from the solution of the 
two-particle problem within a single well. 



Integrating out the single-particle states \l,u) 



"La 



10) 



perturbatively in Ml/vi, we obtain the effective Hamilto- 
nian H M , see Eq. J2J. The last term in Eq. accounts 
for the coupling between molecules \m n ) and the atomic 
states \bi). The operator 6163 — 6264 is the only second- 
order polynomial in bi transforming in the same repre- 
sentation B2 as the d-wave molecule |mn). The coupling 
9 is 



g = -h 2 n 2 
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(6) 



which reduces to g ~ Ah Hi £/o2/(e2 — eo) for weak in- 
teractions C/02 <C eo- We have dropped terms ~ 6? as 
we assume that double occupation of the bosonic sites 
is strongly suppressed by the on-site repulsion U . In 
addition to the Hamiltonian in Eq. J5J), we obtain a 
laser induced hopping term ^}Zi n Jid\ n d lD with Ji = 
— h 2 fl 2 \wi\ 2 /i'i. The sign of this laser induced hopping 
depends on the detuning. In principle it is possible to 
cancel these terms via interference by an additional far- 
detuned Raman transition. Furthermore, we have ab- 
sorbed a shift in the energy of the molecular energy into 
a redefinition of lu. Tuning the Rabi frequency f2 sets the 
energy scale of the coupling g, while the detuning v is 
controlled by the frequency lu of the Raman transition; 
this allows the system to be tuned through a resonance. 

The zero temperature phase diagram of the two- 
dimensional setup shown in Fig. with the Hamilto- 
nian H = H Bli + H M has not yet been studied. However, 
by considering appropriate limits we suggest there is po- 
tential for interesting physics in the intervening regime. 
We let q be the average number of atoms per unit cell 
(i.e. molecular states counted twice), and consider the 
filling q = 1/2. For large positive detuning v 3> g with 
J » g 2 ' /v, the system reduces to the conventional Bose- 
Hubbard model, see Eq. (JJJ, and the bosons establish a 
superfluid phase due to the incommensurate filling. In 
the opposite limit of large negative detuning (\v\ J,g 
and v < 0), all bosons are paired into molecules, and 
we can think in terms of an effective molecular Bose- 
Hubbard model at 1/4 filling. For J <C g, perturbation 
theory in g/\v\ generates a nearest- and next- nearest- 
neighbor repulsion between molecules U M ~ 5 4 /|^ 3 |- Fur- 
thermore, a molecular ring exchange term is generated at 
the same order. The ground state of the resulting model 
is not known, but is likely to be a molecular charge den- 
sity wave with (m^m^) larger on every other row and 
column. There is no difference at the level of symmetry 
between this state and a "plaquette" valence bond solid 
(VBS) of the atomic bosons. We therefore suggest that 
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this system is a candidate for a deconfined quantum crit- 
ical point as the system is tuned between the VBS 
and the bosonic superfluid. 

The above design of a ring exchange interaction for 
bosonic systems is a building block that can be applied 
to many different setups and lattice structures [21|. Of 
particular interest are models exhibiting a local gauge in- 
variance that may allow for the realization of deconfined 
insulators. 

In the following, we focus on the three-dimensional 
setup shown in Fig. ^p, that can likely access such a 
U(l) deconfined state. The lattice structure is described 
by a cubic lattice with the bosonic states on the links. 
Each cubic face forms a plaquette involving four link 
bosons created by the operator b|-, and the molecules 
are again placed at the plaquette centers, i.e. in the cen- 
ter of each cubic face. The lattice of bosonic states can 
also be viewed as the lattice of corner-sharing octahedra 
with their centers at the cubic sites. With this definition 
of the bosonic and molecular sites, we can again derive 
the Hamiltonian H = H BH + H M following the procedure 
discussed above for the two-dimensional setup [2l|. We 
expect that the phase diagram of this system is domi- 
nated by a superfluid phase for large hopping, while in 
the limit of small hopping more exotic states may result. 
To better understand the possibilities, we consider the 
limit of vanishing hopping J = 0. For large detuning 
| v\ ^> <?, the molecules can be integrated out to obtain 
the cubic ring exchange model of Ref. |(| , with an addi- 
tional interaction term, see Eq. 1(3} . (Note, that a unitary 
transformation allows one to change the sign of the ring 
exchange term.) This model is a U(l) lattice gauge the- 
ory and, at least over some region of the parameter space, 
it can enter its Coulomb phase, i.e. a U(l) deconfined 
insulator It has been shown that this state is sta- 
ble to all perturbations, including those that break the 
gauge invariance such as a boson hopping J. Then, in 
the presence of a small but finite hoppping J, the gauge 
structure goes from an explicit, microscopic property to 
an emergent one present only in the low-energy theory. 

Remarkably, the three-dimensional setup shown in 
Fig.Q]with the Hamiltonian H — H BH +H M and quenched 
hopping J = is even a U(l) lattice gauge theory for 
arbitrary gjv. This can be seen by considering a cu- 
bic site i, letting L(i) be the 6 cubic links and P(i) 
the 12 plaquettes containing i. Then the local object 
G % = Ey eL(i) b lj b ij+T,nepty m hm n is a conserved U(l) 
"gauge charge." A straightforward "spin wave" treat- 
ment allows one to write down a low-energy theory of 
liquid phases in ring exchange models [6J, |14|, and we 
obtain here an artificial photon mode with two polariza- 
tions. This analysis demonstrates that this model is likely 
to support a U(l) deconfined insulator, which is likely to 
be continuously connected to the U(l) deconfined phase 
discussed above and in Ref. |6j for the large- v limit. We 
also note that H = H BH + H M should be amenable to 



quantum Monte Carlo simulation when J = 0, as a sim- 
ple unitary transformation renders all matrix elements of 
e --rff nonnegative in the number basis. 

Finally, we remark that the ci-wave molecular state 
couples to a special structure of the surrounding bosons. 
In addition to the applications discussed above for the de- 
sign of strongly correlated models, it should be possible 
to use the coupling of this state as a probe of this struc- 
ture within phases of the single particle system. This 
may provide a powerful tool for the detection of uncon- 
ventional order in strongly correlated systems. 

We are grateful to A. Muramatsu for helpful discus- 
sions. This research is supported by the Austrian Science 
Foundation and the Institute of Quantum Information 
(H.P.B. and P.Z. ), and NSF Grant Nos. DMR-0210790 
and PHY-9907949 (M.P.A.F. and M.H. ). 



[1] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[2] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, and 

I. Bloch, Nature 415, 39 (2002). 
[3] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and 

T. Esslinger, Phys. Rev. Lett. 92, 130403 (2004). 
[4] C. D. Fertig, K. M. O'Hara, J. H. Huckans, S. L. Rol- 

ston, W. D. Phillips, and J. V. Porto, cond-mat/0410491 

(2004). 

[5] D. Jaksch and P. Zoller, Ann. Phys. 315, 52 (2004). 

[6] M. Hermele, M. P. A. Fisher, and L. Balents, Phys. Rev. 

B 69, 064404 (2004). 
[7] R. Moessner, S. L. Sondhi, Phys. Rev. B 68, 184512 

(2003). 

[8] J. B. Kogut, Rev. Mod. Phys. 51, 659 (1979). 
[9] R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 
(2001). 

[10] L. Balents, M. P. A. Fisher, and S. M. Girvin, Phys. Rev. 

B 65, 224412 (2002). 
[11] O. I. Motrunich and T. Senthil, Phys. Rev. Lett. 89, 

277004 (2002). 

[12] D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 

2376 (1988). 
[13] T. Senthil, cond-mat/0411275 (2004). 
[14] A. Paramekanti, L. Balents, and M. P. A. Fisher, Phys. 

Rev. B 66, 054526 (2002). 
[15] A. W. Sandvik, S. Daul, R. R. P. Singh, and D. J. 

Scalapino, Phys. Rev. Lett. 89, 247201 (2002). 
[16] T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and 

M. P. A. Fisher, Science 303, 1490 (2004). 
[17] O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. 

H&nsch, and I. Bloch, Phys. Rev. Lett. 91, 010407 (2003). 
[18] L.-M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. 

Lett. 91, 090402 (2003). 
[19] H. Katori, T. Ido, and M. Kuwata-Gonokami, J. Phys. 

Soc. Japan 68, 2479 (1999). 
[20] L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts, 

H. Fehrmann, and M. Lewenstein, Phys. Rev. Lett. 93, 

030601 (2004). 

[21] H. P. Biichler, M. Hermele, S. D. Huber, M. P. A. Fisher, 
and P. Zoller, in progress. 



